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First- and second-order phase transitions in scale-free networks 
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We study first- and second-order phase transitions of ferromagnetic lattice models on scale- 
free networks, with a degree exponent 7. Using the example of the g-state Potts model we derive 
a general self-consistency relation within the frame of the Weiss molecular-field approximation, 
which presumably leads to exact critical singularities. Depending on the value of 7, we have found 
three different regimes of the phase diagram. As a general trend first-order transitions soften with 
decreasing 7 and the critical singularities at the second-order transitions are 7 dependent. 
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Complex networks, which have more complicated con- 
nectivity structure than periodic lattices (PLs) have at- 
tracted considerable interest recently [1,2]. This research 
is motivated by empirical data collected and analyzed in 
different fields. Small- world (SW) networks [3], which 
can be generated from PLs by replacing a fraction p of 
bonds by new random links of arbitrary lengths, are suit- 
able to model neural networks [4] and transportation sys- 
tems [5]. On the other hand, scale- free (SF) networks [6] 
are realized among others in social systems [7], in pro- 
tein interaction networks [8], in the internet [9] and in 
the world-wide web [10]. In a SF network the degree 
distribution Pr>(k), where k is the number of links con- 
nected to a vertex, has an asymptotic power-law decay 
Po{k) ~ fc -7 , thus there is no characteristic scale in- 
volved. In natural and artificial networks the value of 
the degree exponent is usually in the range 2 < 7 < 3 
[11]. 

Cooperative processes such as spread of epidemic dis- 
ease [12], percolation [13], Ising model [14,15], etc. have 
also been studied in the SW and the SF networks. For 
SW networks numerical studies show [16] that any fi- 
nite fraction of new, long-range bonds, p > 0, brings the 
transition into the classical, mean-field (MF) universality 
class. It is understandable since for systems with long- 
range interactions the MF approximation is exact. In the 
SF networks, where links between remote sites exists, too, 
at first thought one could expect also a traditional MF 
critical behavior. In specific problems, however, it turned 
out that it is only true for losely connected networks, 
when the degree exponent 7 is large enough. Otherwise 
the critical singularities of the transition are model in- 
dependent, but nonuniversal; the critical exponents con- 
tinuously depend on the value of the degree exponent. 
In particular, for 2 < 7 < 3, when (k 2 ) is divergent the 
systems are in their ordered phase for any value of the 
control parameter (temperature, percolation probability, 
transition rate, etc.), and the critical properties can be 
investigated in the limit of infinitely strong fluctuations. 

Till now investigations on cooperative processes in the 
SF networks are almost exclusively limited to continuous 



phase transitions. However, in many problems the phase 
transitions on PLs are first order and it seems natural 
to ask what happens with these transitions on the SF 
networks? There is a general tendency that the discon- 
tinuities (e.g., the latent heat) in the pure system are 
reduced due to inhomogeneities, which often change the 
transition into a continuous one. This has been observed 
in the vicinity of free surfaces [17], when there are miss- 
ing bonds, or in the bulk when random [18] or aperiodic 
[19] perturbations are present. 

In the present paper, we investigate this issue on the SF 
networks. In particular we are interested in the combined 
effect of strong connectivity and irregularities, present 
in the SF networks, on the properties of discontinuous 
phase transitions. In the actual calculations we start with 
the ferromagnetic g-state Potts model and solve it in the 
frame of the Weiss molecular-field approximation, which 
represents a lattice version of the MF method. Then we 
generalize this procedure for any lattice model and show 
how the MF equation on the SF networks can be deduced 
from the corresponding one for PLs. The MF equation 
is analyzed by standard methods [20] and the properties 
of the phase transitions, in particular those related to a 
first- to second-order crossover are calculated. Since the 
MF method is expectedly exact for the SF networks our 
results are presumably exact. 

In the following, we consider the g-state ferromagnetic 
Potts model [21] defined by the Hamiltonian: 



H 



K i] S(s l , Sj) + hiS(si) 



(1) 



in terms of Potts spin variables, Sj = 0, 1, . . . , q — 1, at 
site i. The interaction Kij is equal to K > if the bond 
(ij) is occupied and zero, otherwise. As is well known, 
the Potts model contains as special cases the Ising model 
for q = 2 and the bond percolation problem in the limit 
q — ► 1. On regular, d-dimensional lattices in the absence 
of external fields the phase transition of the homogeneous 
model is first order, as in the MF theory, for q > q c (d) 
and continuous for q < q c (d) where q c {2) — 4, g c (3) < 3 
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and q c (d > 4) = 2. 

To find the thermodynamical properties of the model 
we use the MF method, when the problem is transformed 
to a set of independent spins in the presence of effective 
local fields, which are created by the nearest neighbors. 
The partition function is then given as a product of sin- 
gle site contributions, Z ~JJ i z i , and the free energy F 
takes the form, 



k B T 



2rrij — (q — l)mjmj] 



+ Em 



exp ( ^2 Kij-rrij + hi) + q - 1 



(2) 



Here we introduced the local magnetization < nij < 1 

as 



TO; 



(3) 



the value of which follows from the extremal condition 
of the free energy dF/drrii — 0, leading to a set of self- 
consistency (SC) equations: 

exp (V . Kijrrij + hi) - 1 

£ Karm = ]T K iS ^ . (4) 

i i exp\y2j K ij m j + h i) +9-1 

On a PL with coordination number, z, rrii = m and 
hi = h one obtains the relation, 



m = G(zKm + h), G{x) = G P {x) 



e x - 1 
e x + q — 1 



(5) 



which is compatible with a first-order (second-order) 
transition for q > 2 (q < 2). 

For a SF network we consider no correlations (au- 
tocorrelations) between the degrees of connected sites 
and assume that the probability of having a link be- 
tween sites i and j, Pij is proportional to the number 
of links connected to each sites, i.e., <~ kikj. Fur- 
thermore, in the spirit of the MF method we replace 
the interaction, K i j 1 in Eq. (4) by its average value [22], 
Kij — » K (kikj j J^i ki) ■ Now in terms of the average 
order parameter, to = X)i^i m i/X)i^i an d f° r homo- 
geneous field hi — h one obtains from Eq. (4) the SC 
equation for the SF networks: 



TO 



= J dkP D (k)k G{kKm + h)/(k) = G SF {Km, h) 



(6) 



where summation with respect to i is replaced by integra- 
tion over the degree k as (1/N) ^ — ► J dkPjj{k), where 
N is the number of vertices. Note that the SC equations 
for PLs in Eq. (5) and for the SF networks in Eq. (6) 
are in similar form, and the SC function for networks, 
GsF{Km,h) is directly related to that in a PL, G(x). 
This latter transformation, as given in Eq. (6) remains 



the same for any type of lattice model. Therefore, Eq. (6) 
sets a direct connection between the MF solutions in PLs 
and in the SF networks and thus it is a fundamental re- 
lation. 

Next, we turn to analyze the critical behavior of the 
SF networks compatible with the general SC equation 
in Eq. (6). First, we recall that the SC function, G(x) 
is monotonically increasing from to 1 as x varies from 
to oo and the first few terms of its Taylor expansion, 
G(x) = J2 n =i a nX n are essential for the properties of the 
phase transition [20] . For the Potts model the first three 
coefficients are given by a\ = 1/q, ai = (q — 2)/(2q 2 ), 
and a 3 = (q 2 — 6q + 6)/(6g 3 ). For the SF networks, the 
analogous SC function, Gsf(Kiti, h), is generally not an- 
alytical due to singularities caused by integration over 
the degree distribution. For small m (and for small h) it 
has generally a regular part which is a polynome of finite 
degree n where n is the largest natural number smaller 
than 7 — 2: 



" (k n+1 ) 

G r SF (Km,h) = J2 a "- L Wr( Km ') n + a i h > ( ? ) 

71 = 1 ^ ' 

and a singular contribution, which in the small m limit 
is given by 

_ f a s {Km)i- 2 , n + 2< 7 <n + 3 

^SF(^rn) - < B) ]nKm \( Km -)fi+i i 7 = n + 3 



Here 



■/ 

Jo 



C / dxx 1 -' 1 



G(x) - a « 3 



n=l 



(8) 



(9) 



and the constants B and C are positive. 

Having the small to behavior of the SC function for 
the SF network at hand, we can analyze the correspond- 
ing critical behavior. Due to the presence of the singular 
contribution in Eq. (8) the critical behavior of the SF 
network can be different of that in PLs. Generally we 
can define three regions of the degree exponent 7 with 
different types of critical behavior. In the following, we 
are going to describe these regimes. 

7 > j u : Conventional mean-field regime. If the degree 
exponent is larger than an upper critical value 7" the 
critical behavior on the SF network is the same as on a 
PL. This happens when the singular term in Eq. (8) does 
not modify the usual Landau-type analyzis [20]. Here, 
depending on the order of the transition in the PL, there 
are two different possibilities: 

(i) If the transition in the PL is second order, then 
the first two nonvanishing terms of the small to expan- 
sion of GsF(Km,h) should be regular, i.e., h > 2 and 

Sn=2 l a ™l > 0- In this case 7" = + 2, where n2 > 1 
is the smallest integer for which a„ 2 < 0. As an example 
for the (q < 2)-state Potts model (including percolation 
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for which q = 1) the upper critical degree exponent is 
7" = 4, since here 02 < 0. On the other hand for the 
Ising model, where a 2 = and a 3 < the upper critical 
value is 7" = 5; 

(ii) In the second case, when the transition in the 
PL is first order, we are looking for the condition that 
the transition stays first order in the SF network, too. 
This will happen, provided (a) the linear regular term 
of GsF{Km,h) exists and (b) the next-order contribu- 
tion (cither regular or singular) is positive. It is then 
easy to see that the upper critical value of 7" is given 
by the conditions h = 1 and a s (^ u ) = 0. Indeed, for 
strongly connected networks, when a s < the transition 
is softened into a second-order one, the unconventional 
properties of which will be described later. As an exam- 
ple the first-order transition of the (q > 2)-state Potts 
model in PLs, where 02 > and 03 < 0, will turn into a 
continuous one on the SF networks for 7 < 7", where the 
upper critical exponent obtained numerically, is shown in 
Fig. 1 for different values of q. 




q 

FIG. 1. Regions of first- and second-order phase transitions 
for the g-state Potts model on scale-free networks with a de- 
gree exponent 7. In the second-order regime, i.e., below the 
upper critical value 7" the singularities are 7 dependent. 

As a general trend 7" is monotonously decreasing with 
q and approaching the limiting value of 3 as 1 jq for large 
q (see the inset to Fig. 1). This is consistent with our 
expectations; a stronger first-order transition on a PL, 
which has a larger latent heat, can be destroyed only in 
a more connected network, i.e., with a smaller value of 
7- 

Thus we can conclude at this point that for 7 < 7" the 
effect of the connectivity of the SF network is relevant, so 
that the singularities of the thcrmodynamical quantities 
of the system are different from the conventional mean- 
field behavior, which can be observed in PLs, The rele- 
vant perturbation region is still divided into two parts, 
depending on the position of the singularity: wether it is 
at finite or at zero coupling. In the following, we describe 
these regions. 



3 < 7 < 7": Unconventional critical regime. The criti- 
cal behavior in this regime is due to an interplay between 
a regular linear term (which does exists, since 7 > 3) and 
a negative singular next-to-leading term in the expansion 
of GsF{Km,h). As a result the transition is second or- 
der and takes place at a finite coupling, which in the MF 
method is given by K c = (k) / ((k 2 )ai). Due to the 7 de- 
pendence of the singular term the singularity of the order 
parameter is unconventional: 

m(K) ~ (AK)^^- 3 \ 3< 7 <7 tt , (10) 

where AK = K — K c . At the upper critical value of 
the degree exponent 7 = 7" according to the result in 
Eq. (8), there are logarithmic corrections of the form 

/ ak y/^- 3 ) 

m ^-{WAK\) ' <"> 

at least if the transition in the PL is second order. 

If the transition in the PL is first order, then 7 = 7" 
corresponds to a tricritical point in the SF network and 
the tricritical exponents depend on other details of the 
degree distribution, such as the next-to-leading decay ex- 
ponent. 

The behavior of the susceptibility at the transition 
point is calculated from the small h expansion of the 
SC function in Eq. (7). Since the leading contributions 
are regular, the singularity of the susceptibility follows 
the conventional Curie- Weiss law, x{K) ~ Vl^if |, and 
is not modified by the connectivity effect of the SF net- 
work. 

The singularity in the specific heat is directly related 
to that of the order parameter and can be deduced from 
the known relation for the energy density e ~ m 2 valid 
in the MF theory. 

7 < 3: Ordered regime. If the degree exponent of the 
SF network is 7 < 3 (but 7 > 2, in order to ensure a finite 
average degree, (k) < 00), then the singular properties 
of the system are exclusively determined by the leading 
singular term of the SC function in Eq. (8). As a con- 
sequence the system in the SF network is in its ordered 
phase at any finite value of the coupling and singularities 
take place only at zero coupling (or at infinite tempera- 
ture). The order parameter vanishes at K = as 

m(K) ~ ^-2)7(3-7) ) 2< 7 <3, (12) 

whereas at the borderline value, 7 = 3 there is an essen- 
tial singularity: 

m(K) ~ K^expi-l/BK), 7 = 3. (13) 

The susceptibility at K = is generally finite, \ — 
ai/(3 — 7), except for 7 = 3, when it is divergent as 

In a finite network with N vertices the order in the sys- 
tem disappears already at a nonzero coupling, Kc(N), 
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which can be estimated as follows. The typical value 
of the largest degree in the finite network, k maxi is 
obtained from the usual condition for extreme events: 
f~ P D (k)dk ~ 1/N, thus k max ~ N 1 '^-^. In a finite 
system the different moments of k are also finite, and we 
obtain for the finite-size scaling behavior of the second 
moment: (fc 2 ) - k^ x ~ iV(3-7)/(7-i). From thig rcsult 

the size dependent value of the coupling at the transition 
point can be calculated as 

K C {N) ~ (fc 2 )- 1 ~ 7V-( 3 -7)/(7-i) ) 2 < 7 < 3 , (14) 

from which the finite-size scaling behavior of the order 
parameter at the transition point 

m(N) ~ K c {N)^- 2 V^-i) _ JV -( 7 -2)/(7-i) ( 15 ) 

follows. For 7 = 3 the size dependence of the transition 
point is logarithmic: 

Kc{N)~(\nN)~\ 7 = 3, (16) 

which has been observed in Monte Carlo simulations for 
the Ising model [14]. 

To summarize we have studied the properties of first- 
and second-order phase transitions of ferromagnetic lat- 
tice models on scale-free networks. Using the Weiss MF 
approximation we have derived a general SC equation 
that has been analyzed by standard methods. Three re- 
gions of the phase diagram are identified as a function of 



the degree exponent 7. First-order transitions can only 
take place in the first regime, where 7 > 7", and where 
the effect of connectivity of the network is irrelevant. In 
the second regime, for 3 < 7 < 7" the phase transition 
is always continuous and takes place at a finite temper- 
ature. The magnetization critical exponent, however, is 
nonconventional, 7 dependent. In the third regime, for 
7 < 3, the system is in its ordered phase at any tem- 
perature and a singularity can develop only for diverging 
temperature (or vanishing coupling). 

As far as the properties of the critical singularities are 
concerned the MF method presumably gives exact re- 
sults. The location of the transition point is not neces- 
sarily exact. Indeed, our results coincide with others, ob- 
tained by different methods on specific problems [13,15]. 

After this work has been completed we became aware 
of a preprint by Goltsev, Dorogovtsev, and Mendes [23], 
in which some results about the (/-state Potts model on 
the SF networks have been announced. 
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